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Aims. Our goal is to research the upper limits on the extragalactic background light (EBL). 

Methods. The upper limits on the extragalactic background light (EBL), using the Fermi and very high energy (VHE) spectra recently 
observed in TeV blazars, are presented. We use an assumption that the VHE intrinsic photon index cannot be harder than the Fermi 
index measured by the Fermi-LAT. 

Results. Totally, these upper limits on the EBL are compatible with ones given by most of EBL models. However, the models of high 
EBL density are denied by TeV blazars. 
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1. Introduction 

The diffuse extragalactic background light (EBL) consists of the 
sum of the starlight emitted by galaxies through the history of 
the Universe, and includes an important contribution from the 
first stars. Direct measurements of the extragalactic background 
light (EBL) from the infrared (IR) to the ultraviolet (UV) are dif- 
ficult because of the light pollution of bright foreground sources 
(see comprehensively reviewed measurements and implications 
of the cosmic infrared background, Hauser & Dwek 2001). The 
method of galaxy counts is used to estimate the EBL, it just pro- 
vides a lower limit owing to the unknown of unresolved sources 
Various models for the EBL have been published (Salamon & 
Stecker 1998; Malkan & Stecker 1998; Malkan & Stecker 2001 
Stecker et al. 2006; Kneiske et al. 2002; Kneiske et al. 2004 
Primack et al. 2005; Primack et al. 2008; Gilmore et al. 2008 
Gilmore et al. 2009; Franceschini et al. 2008; Razzaque et al 
2009; Finke et al. 2010). These models include different degrees 
of complexity, observational constraints and data inputs. 

Absorption features imprinted on the very high en- 
ergy (VHE) spectra of distant extragalactic objects by 
background light photons provide an indirect approach to 
study the EBL. Assuming an intrinsic gamma-ray spec- 
trum, [Stecker & de Jager( 1993)1 |Stanev & Franceschini(1998)| 
have constrained the EBL from the observed VHE spectra of 
blazars. p\.haronian et al.(2006a)| have also discussed upper lim- 
its on the background light at optical/near-infrared wavelengths 
based on the HESS observation of lES 1101-232. They as- 
sume the intrinsic spectrum to be not harder than r,„, = 1.5 
and put limits on EBL quite close to the lower limits by galaxy 
counts. The detail studies of EBL shapes and blazar VHE spec- 
tra are also given by [Mazin & Raue(2007)| |Schroedter(2005)| 
|Fink e & Razzaque(2009)' where same intrinsic spectrum for all 
blazars is assumed. In fact, blazars have different intrinsic spec- 
tra. The main handicap of this approach to limit EBL is the un- 
certainty about the intrinsic spectrum of VHE. 



To date, 35 AGN sources have been detected at TeV 
energies(E> 1 00 Ge Vfl Their observed VHE spectra have power 
law shapes with the index Fvhe ^ 2, in which distant sources 
have large Vvhe, up to 4 (e.g. Acciari et al. 2009b; Albert et al. 
2007b, 2008b; Aharonian et al. 2006c, 2005a). Many of these 
sources have recently been detected at GeV energies by the 
Fermi Gamma-ray Space Telescope(Abdo et al. 2009; 2010b). 
Ab do et al. (20 09) have extrapolated the Fermi spectrum up to 
10 TeV assuming a single spectral index and taken it as an in- 
trinsic spectrum of VHE ranges. The most break of the observed 
VHE spectra are consistent with the absorption predicted by the 
minimal EBL density model. For a TeV source, the presence 
of a break between the Fermi and VHE energy range might 
be caused by some internal or external factors. The internal 
factors include a break of emitting particle distribution or an 
intrinsic absorption caused by strong optical-infrared radiation 
within the source (Donea & Protheroe 2003). The external fac- 
tors usually refer to the cosmic attenuation effect. Furthermore, 
it is difficult to well predict the intrinsic spectrum from simul- 
taneous multi-wavelength observations because of the complex- 
ity of VHE emission mechanism. In this work, we assume that 
the Fermi spectral index measured by Fermi-LAT is the lower 
limit of intrinsic VHE spectral index for TeV blazars instead 
of single 1.5. In the other words, the photon index from the 
Fermi to VHE energy range is only softened except the pres- 
ence of a new component (Yang & Wang 2010), or monochro- 
matic radiation fields within the source (Aharonian et al 2009a). 
We note that the steeper intrinsic index assumed by us provides 
stronger constraints to the EBL intensity (all Fermi photon in- 
dex of TeV sources larger than the 1 .5 with the exclusion of H 
1426-1-428 Tfer - 1-47). Moreover, taking into account the dif- 
ferences of VHE emission between the different sources, the as- 
sumption is more reasonable than one of > 1.5. Recently 
'Georganopoulos et al. (2010)'and 'Mankuzhiyil et al.(2010) also 
use the extrapolation of the Fermi data as upper limits of in- 
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trinsic TeV spectra. Assuming the VHE intrinsic spectra cor- 
rected by EBL absorption be softer than the Fermi spectra, 
[Prandini et al. (201 0) found that the derived redshifts are larger 
than true ones. It shows that their assumption is reliable. Based 
on this assumption, we analyze the Fermi and VHE spectra of 
TeV blazars and give upper limits on the EBL intensity. 

In § 2 we describe the method of calculating y - y absorp- 
tion optical depth, Tyy(e), and the EBL intensity, developed by 



the energy e' = lin^c'^ /(Ey{l + z)) ~ 2m^c'^/Ey where the pair- 
production cross section reaches the largest value, and give an 
upper limit on the EBL number density, n(e, z). Using the Dirac 
delta-function, we approximately write n(e, z) as 



|Schroedter(2005)| and [Fmke & Razzaq ue(2009)i In § 3 we apply 
the method to the TeV blazars with VHE and Fermi spectra and 
discuss the limits of these sources on the EBL. 



2. THE METHOD 

The VHE absorption of the EBL is caused by the pair production 
of photon-photon collision. The observed VHE flux is given by 



n{£,z ~0)~ 6'n(e',z ^ 0)6(e-e'). 
Integrating Eq.Q, we obtain: 

2H[)T'yy''Ey 



n{e',z ~ 0) 



cZe^r^rrP-C^ifiil) 



(6) 



(7) 



where (f{2) ^ 1.787. The error of the EBL number density is 
given by 



o-(n) = 



2Ho(r(T^^')Ey 



(8) 



Us(Ey) = e-'^^y^fi„,{Ey) , 



(1) 



where T(Ey) is the optical depth, Ey is the observed y-ray photon 
energy, and fi,„iEy) is the intrinsic flux. 

For a VHE source at redshift Ze, the optical depth of its Ey 
energy photon caused by the EBL is given by 



T(Ey,z<,) = cnrJ- 



Wz^ r ^ 

dz J^»kL 



de-e ^n(e,z)^[io(e)](2) 



where n(e,z) is the photon number density of the EBL 
with energy e at redshift z, is the classical electron ra- 
dius. So = eEylm^c^, iplsoie)] is a function given by 
[Gould & Schreder( 1967)1 and ^ is the diff'erential time of red- 
shift given by 



dt/dz 



\(l + zf(l + n„.z) - z(z + 2)Qa1"'^' . 



(3) 



|Abdo et al. (2009)] have found that the intrinsic spectra of 
many TeV sources can be described by a single power-law across 
the Fermi and VHE energy ranges. In fact, the observed Fermi 
and VHE spectral indices are different duo to the EBL absorp- 
tion. Their diff'erence AF increases with redshift. For example, 
M 87 and Cen A with low redshidts have AF x 0, while blazars 
with redshifts greater than 0.1 show AF > 1.5. We assume that 
the observed Fermi spectral index is the lower limit of intrinsic 
VHE spectral index, i.e., T''"J"y^^ ^ ^per- For some objects, mul- 
tiple VHE spectra have been observed in different flux states. 
We adopt the VHE spectra of their low flux states to constrain 
the EBL. 

Based on r"™^,^^, an upper Umit on the optical depth 
T(Ey, Ze) is given by |Finke & Razzaque(2009)| 

7^ {Ey^Ze) — '^{Ey ifiifi, Ze^ + (J^obs ^ int VHE^^^^^y ^ ^V^^'^'^^^ 

and its standard error t'""^ is given by 

cr(T;7) = CT{UAEy))lfobAEy), (5) 



cZe''^rlm^c'^^p{2) 
Finally the EBL intensity is given by 

vIAz) ^ ^ Me,z) . (9) 
At: 

3. Results and Discussion 

The EBL has two spectral humps with different origins. The blue 
hump at UV-Optical-NIR (near-infrared) wavelengths comes 
from stars. The red hump at MIR (mid-infrared) and FIR (far- 
infrared) wavelengths is from the absorption and re-emission of 
starlight by the interstellar medium. Therefore, the EBL includes 
the important information of star formation and evolution. 

The TeV blazars used to constrain the EBL are listed in 
Table 1, where the spectra at low-flux state are used to the ut- 
most, since the 1 1 months averaged Fermi spectra are unlikely 
to correspond to the high state. Through calculation, we find that 
four blazars, 3C 66A, 0716-1-714, 3C 279, and PG 1553+113, 
give stronger constraint on the EBL density. Since other TeV 
blazars are consistent with all listed EBL models within the er- 
ror range, we do not give their constraint. The EBL upper limits 
given by the spectra of four TeV blazars are shown in Fig. [T] 
The curves of several EBL models are also plotted in Fig. [1] 
Gilmore et a l.(2009) |Stecker et al.(2006)| 



Kneiske et al.(2004) 



Finkeet al (2010) 



ranceschini et al.(2008) In this work the 



where T(Ey^,„i„,Ze) at the lowest energy -Ey,,,,,-,, of 
VHE observations is estimated by the EBL model of 
[Franceschini et al.(2008)l 

Now we use ■f™'^(Ey,Ze) to estimate an upper limit 
on the EBL number density. Following ' S chroedter(2005 )[ 



Finke & Razzaque( 2009)| we take dtldz ~ //q ' for TeV sources 
due to low redshift, where Hq - 70 km s ' Mpc ^ We assume 
the monochromatic absorption of VHE photon Ey by the EBL at 



calculated limits on the EBL at UV-Optical-NIR wavelengths 
are strong. For comparison, we also list the lower limits of the 
EBL from source counts [Madau & Pozzetti(2000)| shown by the 
upward triangles. These calculated limits are inconsistent with 
the fast evolution model given by [Stecker et al.(2006)| in NIR- 
Optical wavelengths, but are still compatible with their baseline 
model. For the fast evolution model of St ecker et al.(2006)| the 
extinction of UV photons by the interstellar gas in galaxies is 
not considered, the UV and Optical-NIR photon density might 
be over-estimated. In fact, the observed gamma-ray hard spectra 
of H 2356-309 (z = 0.165) and lES 1101-232 (z = 0.186) by 
[Aharonian et al.(20 06a) suggest that an upper limit to the EBL 
at optical- NIR wavelengths is very close to the lower limit given 
by the integrated light of resolved galaxies. This implies that 
the EBL is more transparent to high energy y-rays than previ- 
ously thought and the contribution from sources except starlight 
is less. |Essey & Kusenko(2009) suggest a new interpretation of 
these observations. For distant blazars, the gamma-ray emis- 
sion is dominated by the secondary photons, while for nearby 
blazars the emission is from the primary photons. Therefore, 
they argue that distant AGN would show no significant atten- 
uation due to pair production on the EBL. Also, we should 
note that the redshift of 3C 66A, assumed to be z = 0.444, 
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Fig.l. Upper limits of the EBL given by 3C 66A, 0716+714, 
3C 279 and PG 1553+1 13 under the assumption that the Fermi 
spectral index measured by Fermi-LAT can be used as an lower 
limit of the intrinsic VHE spectral index. Also plotted are sev- 
eral EBL models: the baseline and fast evolution models of 
[Stecker et al.(2 006ydash, dot green curves, respectively), the 
model of [Franceschini et al.(2008Xdash red curve), the fiducial 
model from Gilm ore et al. (2 009) 'dot black curve), the best fit 
model from Knei ske et al.(2 004)vdot-dash blue curve), and the 
Model C of Fink e et al. (201 0)vsolid black curve). We also Hst 
the lower limits of EBL from source counts with upward trian- 
gles (Madau & Pozzetti 2000). 



has large uncertain (Miller et al. 1978). If its redshift is less 
than 0.444, the upper limits obtained will be relaxed. Due to 
the very bright nucleus, the redshift of 0716+714 is still un- 
certain. Stickel et al. (1993) repeated spectroscopic observations 
without obtaining the redshift, but they found that two neigh- 
boring galaxies have quite similar redshifts of 0.26. However, 
its host galaxy detection gave the redshift of z = 0.31 (Nilsson 
et al. 2 008). The redshift of P G 1553+113 also has large un- 
certain.|Sbarufatti et al. (2005)| derived its lower limit of z>0.78. 
[Sbarufatti et al. (2006) used the spectra of ESO VLT to give a 
limit of z>0.09.|Abdo et al. (20 10a)| constrained z < 0.75 com- 



bining Fermi and VHE gamma-ray data. Danforth et al.(2010) 
constrained its redshift to be z ~ 0.4 - 0.6 by Hubble Space 
Telescope. In this work, we adopt its redshift as z=0.78. While 
3C 279 has well-known redshift. Totally, in spite of large redshift 
uncertainty for some objects, the models of high EBL density are 
denied by 3C 279. 

In low redshift, we assume the EBL do not evolve with red- 
shift. However, the EBL is progressively generated by galaxies 
and active nuclei (AGN) during most of the Hubble time, partic- 
ularly below z = 1 . The evolution of their photon number density 
is a very complex function of time and frequency. 

The main problem of TeV blazars observations limiting on 
the EBL is the uncertainty of VHE intrinsic spectra. In the 
early works, all TeV blazars are assumed to have a lower limit 
of Tint - 1.5 based on some theories of particles accelera- 
tion and emission (e.g. Aharonian et al 2006a). The same as- 
sumption was also adopted by S chroedter(2005 )| with the anal- 
ysis of observed correlation between their Fvhe and redshift. 



Fig. 2. Difference, AF, between the measured VHE and Fermi 
photon indices (or conventional limit 1.5) as a function of the 
redshift. Red inverted triangles denote the rvHE-^Fer, and black 
triangles denote the Tvhe - 1 .5. The red dash line show the fitting 

of r VHE - f^Fer- 



in the energy ranges of 0.2 TeV < E < 2 TeV, and found that 
Ar(z) - FvHE - r„« is a linear function of the redshift z in the 
range of 0.05 - 0.4. In this work, we assume the Fermi spectral 
index measured by Fermi-LAT to be the lower limit of intrinsic 
VHE spectral index. In Fig.|2] we compare the redshift variation 
of rvHE - Tfe,- and Tvhe - 1 .5 in the range of 0.03 < z < 0.4 ex- 
cluding 3C 66A (z=0.444), 3C 279 (z=0.536) and PG 1553+1 13 
(z=0.78). It is shown that they have similar correlation with red- 
shift, but the parameters fitting TvHE-f'Fer are more similar with 
C and D of the baseline model given by |Stecker & Scully (2006)| 
and Stecker & Scully (2010) In fact, the baseline model gives 
lower EBL density than the fast evolution model does, high 



Stecker & Scully (2006) have derived a simple analytic approx- 



imation of the EBL absorption on the spectra of TeV blazars 



EBL density models are disfavored by observations, such as 
Aha ronian et al.(2006a)[ |Georganopoulos et al. (20l0)| Most of 
the FvHE - Tfer are less than the Tvhe -1-5, and smaller AF will 
provide stronger constraints to EBL. If using r,„, > 1 .5 to limit 
the EBL is reasonable, F,-,,, > Ff^r will be more feasible. 3C 66A 
and 0716+714 obviously deviates the correlation shown in Fig. 
|2] it implies that the assumed redshift might be uncorrect. 

Since the simultaneous data of Fermi and VHE are less avail- 
able nowadays, we only use non-simultaneous spectra. We also 
note that the upper limits of EBL density are great depended on 
the VHE photon index (see the equation In fact, no signifi- 
cant spectral variability is observed in VHE bands. For example, 
PG 1553+1 13 has many times VHE observations (Aharonian et 
al. 2006b, Aharonian et al. 2008, Albert et al. 2007f, Albert et 
al. 2009), its photon index is very similar as observed by HESS 
and MAGIC (Abdo et al. 2010a). For PKS 2155-304, no sig- 
nificant spectral variability appears despite flux variation with 
a factor of two (A haronian et al. 2009). For lES 1218+304, 
[Acciari et al.(2010 ) found that the VHE spectral shape has no 
change between flare and quiescent state. Therefore our calcu- 
lated results are meaningful when VHE spectra have no obvious 
variation. 

We use the assumption of monochromatic absorption to 
calculate the EBL intensity at specified wavelengths. The re- 
sult will be larger than the actual case due to the EBL 
photons near specified wavelength also contributing the ab- 
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Table 1. TeV Blazar Sample 



Blazar 


Redshift 


Fermi photon Index (Tf„) 


VHE photon IndexiTvHE) 


E,„i„ [TeV] 


^max [TcV] 


Reference 


3C 66A 


0.444 


1.93 ±0.02 


4.1 ±0.4 


0.23 


0.47 


1 


S5 0716+714 


0.31 


2.15 ±0.03 


3.45 ±0.54 


0.18 


0.68 


2 


lES 0806+524 


0.138 


2.1 ±0.1 


3.6 ± 1.0 


0.31 


0.63 


3 


lES 1011+496 


0.212 


1.93 ±0.04 


4.0 ±0.5 


0.15 


0.59 


4 


Mai-k421 


0.031 


1.81 ±0.02 


2.2 ± 0.08 


0.13 


2.86 


5 


Mai-k 180 


0.046 


1.86 ±0.11 


3.3 ± 0.7 


0.18 


1.32 


6 


lES 1218+304 


0.182 


1.7 ±0.08 


3.08 ± 0.34 


0.19 


1.44 


7 


W Comae 


0.102 


2.06 ± 0.04 


3.81 ±0.35 


0.27 


1.15 


8 


3C 279 


0.536 


2.32 ± 0.02 


4.11 ±0.68 


0.08 


0.47 


9 


H 1426+428 


0.129 


1.49 ±0.18 


3.5 ±0.35 


0.82 


5.66 


10 


PG 1553+113 


0.78 


1.66 ±0.03 


4.0 ±0.6 


0.21 


0.50 


11 


Mark 501 


0.034 


1.85 ±0.04 


2.45 ± 0.07 


0.14 


4.58 


12 


lES 1959+650 


0.048 


2.1 ±0.05 


2.58 ±0.18 


0.19 


2.40 


13 


PKS 2005-489 


0.071 


1.9 ±0.06 


4.0 ± 0.4 


0.23 


2.27 


14 


PKS 2155-304 


0.117 


1.91 ±0.02 


3.32 ± 0.06 


0.23 


2.27 


15 


BL Lacertae 


0.069 


2.38 ± 0.04 


3.6 ±0.5 


0.16 


0.70 


16 


lES 2344+514 


0.044 


1.57 ±0.12 


2.95 ±0.12 


0.19 


4.00 


17 



The Fermi data c ome from the |Abdo et al. (2010b)| V HE data refer to: (l^Acciari et al.(2009a)| (2) [Anderhub et al.(2009)"l (3) 
lAcciari et al.( 2009b)| ( 4)|Albert et al.(2007b)| (5)|Albert et al.(2007c); (6)IAlbert et al.(2006)| (7)|Accia ri et al.(2009c)[ (8)|Acciari et al. (2008)1 
(9) Albert eta L(2058)l (10) Aharonian et al.(2003) , (11) Aharonian et al.(2006b) , (12) Albert et al.(2Q07d)| (13) [Tagliaferri et al.(2008)[ (14) 
[Aharonian et al.(2005a)| (15) ,Aharonian et al.(2005b), (16) ,Albert et al.(2007a), (17) ,Albert et al.(2007e)| . 



sorption. But, this method proposed by | Schroedter(2005 )| 
Finke & Razzaque(2009)| does not assume the EBL spectrum 
and has the advantage compared to other methods because the 
EBL spectrum is not easily known. Future observations of Fermi 
and VHE spectra for blazars will provide strong constraints on 
the EBL. 
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